arXiv: 1509.0517Ivl [astro-ph.GA] 17 Sep 2015 


Astronomy & Astrophysics manuscript no. guvics lfall rev © ESO 2015 

September 18, 2015 


The GALEX Ultraviolet Virgo Cluster Survey (GUViCS). VI: The UV 
luminosity function of the Virgo cluster and its surrounding 

regions^ 

A. Boselli', S. Boissier\ E. Voyer^ L. Ferrarese^, G. Consolandi^, L. Cortese"^, P. Cote^, J.C. Cuillandre^, G. Gavazzi^, 

S. Gwyn^, S. Heinis^, O. Ilbert\ L. MacArthur^, Y. Roehlly^ 

* Aix Marseille Universite, CNRS, LAM (Laboratoire d’Astrophysique de Marseille), UMR 7326, F-13388, Marseille, France e-mail: 
alessandro.boselli@lam.fr, samuel.boissier@lam.fr, elysse.voyer@gmail.com, olivier.ilbert@lam.fr, 
yannick.roehlly@lam.fr 

^ NRC Herzberg Astronomy and Astrophysics, 5071 West Saanich Road, Victoria, BC, V9E 2E7, Canada e-mail: 

laura.ferrarese@nrc-cnrc.gr.ca, patrick.cote@nrc-cnrc.gr.ca, Stephen.gwyn@nrc-cnrc.gr.ca 
^ Universita di Milano-Bicocca, piazza della scienza 3, 20100, Milano, Italy e-mail: guido.consolandi@mib.infn.it, 
giuseppe.gavazzi@mib.infn.it 

Centre for Astrophysics Supercomputing, Swinburne University of Technology, Mail H29-PO Box 218, Hawthorn VIC 3122, 
Australia e-mail: lcortese@swin.edu.au 

^ CEA/IRFU/SAP, Laboratoire AIM Paris-Saclay, CNRS/INSU, Universit Paris Diderot, Observatoire de Paris, PSL Research 
University, F-91191 Gif-sur-Yvette Cedex, France e-mail: jc. cuillandre@cea. fr 
® Department of Astronomy, University of Maryland, Stadium Drive, College Park, MD 20742-2421, USA e-mail: 
sebastienheinis@gmal.com 

’ Department of Astrophysical Sciences, Princeton University, Princeton, N1 08544, USA e-mail: lauren@astro. princeton. edu 


ABSTRACT 

We use the GALEX data of the GUViCS survey to construct the NUV luminosity function of the Virgo cluster over ~ 300 deg.^, an 
area covering the cluster and its surrounding regions up to ~ 1.8 virial radii. The NUV luminosity function is also determined for 
galaxies of different morphological type and NUV - i colour, and for the different substructures within the cluster. These luminosity 
functions are robust vs. statistical corrections since based on a sample of 833 galaxies mainly identified as cluster members with 
spectroscopic redshift (808) or high-quality optical scaling relations (10). We fit these luminosity functions with a Schechter function, 
and compare the fitted parameters with those determined for other nearby clusters and for the field. The faint end slope of the Virgo 
NUV luminosity function (a = -1.19), here sampled down to ~ NUV = -11.5 mag, is significantly flatter than the one measured in 
other nearby clusters and similar to the field one. Similarly M* = -17.56 is one-to-two magnitudes fainter than measured in Coma, 
A1367, the Shapley supercluster, and the field. These differences seem due to the quite uncertain statistical corrections and the small 
range in absolute magnitude sampled in these clusters. We do not observe strong systematic differences in the overall NUV luminosity 
function of the core of the cluster with respect to that of its periphery. We notice, however, that the relative contribution of red-to-blue 
galaxies at the faint end is inverted, with red quiescent objects dominating the core of the cluster and star forming galaxies dominating 
beyond one virial radius. These observational evidences are discussed in the framework of galaxy evolution in dense environments. 

Key words. Galaxies: clusters: general ; Galaxies: clusters: individual: Virgo; Galaxies: evolution; Galaxies: interactions; Galaxies: 
ISM; Galaxies; star formation; 


1. Introduction 

The environment plays a fundamental role in shaping galaxy 
evolution. Since the seminal work of A. Dressier in the eighties 
it became evident that galaxies in rich environments are sys¬ 
tematically different than those located in the field. They are 
generally quiescent ellipticals and lenticulars, in contrast to the 
dominant population of star forming disc galaxies along the 
cosmic web (Dressier 1980; Binggeli et al. 1988; Whitmore 
et al. 1993; Dressier et al. 1997). At the same time it became 
clear that the physical properties of the star forming systems in¬ 
habiting rich environments such as clusters and groups are also 
systematically different than those of their isolated analogues. 


* Table A. 1 is available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr(130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ 


with a reduced gas content and activity of star formation (e.g. 
Boselli & Gavazzi 2006). 

Several physical processes have been proposed to explain 
the origin of these differences. They can be divided in two 
main families, those related to the gravitational interactions 
between galaxies or with the potential well of the cluster (tidal 
interactions - Merritt 1983; Byrd & Valtonen 1990, harassment 
- Moore et al. 1998), and those exerted by the hot and dense 
intracluster medium on galaxies moving at high velocity within 
the overdense region (ram pressure stripping - Gunn & Gott 
1972, viscous stripping - Nulsen 1982, thermal evaporation - 
Cowie & Songaila 1977, starvation - Larson et al. 1980). These 
processes might start to become active well before galaxies 
entered rich clusters because these large dynamically bounded 
structures have been formed through the accretion of smaller 
groups of galaxies (pre-processing; Dressier 2004). 




2 


Boselli et al.: GUViCS: the UV luminosity function of the Virgo cluster 


The identification of the dominant environmental process 
affecting galaxy evolution can be done through the statistical 
comparison of the physical properties of galaxies in a wide 
range of densities, from the core of rich clusters, compact and 
loose groups, down to the field, cutting the studied samples in 
slice of redshift to trace the time-evolution of the effects. Indeed, 
these processes act on different timescales and are characterised 
by efficiencies that depend on the density of galaxies and of 
the intergalactic medium, the velocity dispersion within the 
structure, and the mass of galaxies. The identification of the 
perturbing process can also be done through the comparison of 
the radial and kinematic properties of the perturbed galaxies 
with the predictions of tuned models and simulations. 

Thanks to its proximity (~ 17 Mpc, Gavazzi et al. 1999; Mei 
et al. 2007) and to its unrelaxed structure indicating its still un¬ 
dergoing formation, the Virgo cluster has become the ideal target 
in the local universe for the study of the effects of the environ¬ 
ment on galaxy evolution. The advent of large panoramic de¬ 
tectors made possible the observations of the whole Virgo clus¬ 
ter region in the optical bands by the NGVS survey (Ferrarese 
et al. 2012), in the 21 cm HI line by ALFALFA (Giovanelli et 
al. 2005), and in the infrared by Herschel (The HeViCS survey, 
Davies et al. 2010). We have recently completed a UV survey 
of the cluster over ~ 300 deg.^ using the GALEX space obser¬ 
vatory (the GUViCS survey; Boselli et al. 2011). UV data are 
of particular interest to trace the content and distribution of the 
young stellar populations in late-type galaxies (e.g. Kennicutt 
1998; Boselli et al. 2009), while the one of very evolved stars 
in quiescent systems (O’Connel 1999; Boselli et al. 2005a). For 
these reasons the UV data are critical to trace short-term varia¬ 
tions in the star formation activity of late-type systems as those 
encountered in galaxies freshly infalling in rich clusters (e.g. 
Boselli et al. 2006, 2008a,b, 2014a,b; Hughes & Cortese 2009; 
Cortese & Hughes 2009). The UV data extracted from this sur¬ 
vey have been recently published in Voyer et al. (2014), and the 
study of the physical properties of Virgo cluster galaxies derived 
from a multifrequency analysis in Boselli et al. (2014a). In two 
other companion papers we presented the study of the extrapla- 
nar star formation activity of a few gas stripped galaxies in the 
cluster (Boissier et al. 2012), and of the diffuse Galactic emis¬ 
sion of the UV scattered light in the Virgo direction (Boissier 
et al. 2015). In the present work we derive and study the UV 
luminosity function of the Virgo cluster and its surrounding re¬ 
gions. This work is thus an extension of the study presented in 
Boselli et al. (2011), where the UV luminosity function was de¬ 
rived only for the central 12 deg.^. The luminosity function is a 
statistical tool widely used to study the mean properties of large 
samples of galaxies. Often used in other wavelengths domains, it 
has been derived at rest frame in the GALEX UV bands only for 
the general field (Wyder et al. 2005) and for a few other nearby 
clusters (Coma, Cortese et al. 2008; Hammer et al. 2012; A1367, 
Cortese et al. 2005; Shapley, Haines et al. 2011). Compared to 
these previous works, the UV luminosity function derived from 
the GUViCS survey has the major advantage of including a sig¬ 
nificant larger number of objects with direct redshift measure¬ 
ments, of sampling a much wider range in galaxy density, from 
the core of the cluster out to the periphery (~ 2 /?v,,>), and of in¬ 
cluding the dwarf galaxy population {NUV < -11.5 mag). Dwarf 
systems are of particular interest in the study of the effects of the 
environment on galaxy evolution since they are fragile objects 
easily perturbed just because of their shallow potential well (e.g. 
Boselli & Gavazzi 2014). 

The paper is structured as follows: in section 2 we describe the 


dataset used in this analysis, while in section 3 we explain the 
different techniques used to identify Virgo cluster members. The 
UV luminosity function is derived in section 4 and analysed in 
the framework of galaxy evolution in rich environments in sec¬ 
tion 5. 

2. The data 

The GUViCS survey has covered the Virgo cluster in the region 
between 12h < R.A < 13h and 0" < dec < 20°. In this region the 
sky coverage is 94 % (Boselli et al. 2014a) at the depth of the 
AIS, where the completeness is 20 magnitude, and 65 % at the 
depth of the MIS (21.5 magfl Unfortunately, because of a fail¬ 
ure of the EUV detector, this survey has been completed only in 
the NUV. The analysis presented in this work will thus be limited 
to this band. The number of NUV detected sources in the cov¬ 
ered region is more than 1.2 millions (Voyer et al. 2014). These 
are Milky Way stars, nearby galaxies and background sources. 
To construct the NUV luminosity function of the Virgo cluster 
and of its surrounding regions we have to count the number of 
galaxies in a given volume down to a limiting magnitude. Eor 
consistency with our previous works we consider as Virgo mem¬ 
bers those galaxies with a recessional velocity vel < 3500 km 
s“'. This limit has been chosen to include all the objects belong¬ 
ing to the different cluster substructures and their surrounding re¬ 
gions (e.g. Boselli et al. 2014a). As limiting magnitude we adopt 
NUV-20 mag. Down to this magnitude, the catalogue is 100% 
complete over the ~ 300 deg.^ covered by the survey. 

Voyer et al. (2014) provide integrated magnitudes for ex¬ 
tended sources (UV-VES catalogue) included in the widely 
known catalogues of nearby galaxies (NGC, UGC, IC, VCC, 
...), for a total of 1770 galaxies. These magnitudes have been 
extracted with the EUNTOOLS analysis package on DS9 using 
manually defined apertures adjusted to fit the full UV profile of 
each traget galaxy. The same procedures have been adopted for 
all the sources identified as Virgo cluster members using spec¬ 
troscopic redshifts. These procedures are necessary for extended 
sources where the automatic GALEX pipeline often does not 
recognise the extended nature of the galaxies from the point¬ 
like structure of HII regions, galaxy nuclei or complex high sur¬ 
face brightness features present in star forming systems. For the 
other sources, NUV fluxes have been taken from the catalogue 
of point-like objects (UV-VPS) of Voyer et al. (2014). Down to 
the limiting magnitude of NUV-20 mag this catalogue includes 
33026 objects. These are generally faint objects, where the NUV 
magnitudes extracted from the pipeline are accurate enough for 
the purpose of this work. Indeed, as clearly shown in Boselli 
et al. (2011) (Fig. 7), the difference between extended magni¬ 
tudes and those extracted from the standard GALEX pipeline op¬ 
timised for point-like sources is <0.3 mag whenever NUV >17 
mag. UV magnitudes are corrected for Galactic attenuation us¬ 
ing the Schlegel et al. (1998) map combined with the Fitzpatrick 
& Massa (2007) extinction curve. Given the high Galactic lati¬ 
tude of the Virgo cluster, these corrections are low (A(NUV) < 
0.4 mag; Boissier et al. 2015). 

For all galaxies at the distance of Virgo (vel < 3500 km s ') 
we also extracted 22 //m WISE fluxes necessary for the dust in¬ 
ternal attenuation corrections. These fluxes have been extracted 
from the original WISE fits images with the same procedures 
used for the UV data, consistently with Boselli et al. (2014a). 

' The coverage is less than 100 % because of the presence of a few 
bright young stars which have been expressly avoided during the survey 
since saturating the UV detector 
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The dust extinction correction is done following the prescription 
of Hao et al. (2011). NUV corrected data are used to estimate 
the star formation rate in star forming galaxies using again the 
prescription of Hao et al. (2011). 

Distances are necessary to calculate absolute magnitudes. 
For consistency with our previous works, distances are deter¬ 
mined assuming the mean distance of the different cluster sub¬ 
groups (see Table 3 in Boselli et al. 2014a). Galaxies at the pe¬ 
riphery of the cluster are all assumed at 17 Mpc, the distance of 
the main body of Virgo. 

3. The identification of Virgo members 

The construction of the NUV luminosity function needs first the 
identification, among the UV detections with NUV < 20 mag, of 
all galaxies at the distance of Virgo. This is done using spectro¬ 
scopic redshift whenever available (see sect. 3.1), while photo¬ 
metric relations (sect. 3.2) or statistical considerations otherwise 
(sect. 3.3). 

3.1. Identification based on spectroscopic data 

Voyer et al. (2014) has cross-matched all the UV detections with 
the spectroscopic catalogues available on the net (NEE0> SDSS 
DRSo). For these objects the spectroscopic information can be 
easily used to identify Virgo members {vel < 3500 km s“'). Out 
of the 1770 galaxies listed in the extended source catalogue, 688 
galaxies match these criteria. There are, however, 34 extended 
galaxies with NUV < 20 mag but without any redshift informa¬ 
tion. Among these, we identify the possible cluster members us¬ 
ing the same procedures adopted for the other galaxies included 
in the UV point source catalogue (see below). 

As previously mentioned, down to NUV < 20 mag the 
majority of the point-like detected objects are not Virgo clus¬ 
ter members but rather Galactic stars and background sources. 
Given the poor angular resolution of the GALEX images (~ 4-5 
arcsec), the separation of point-like stars from extended galaxies 
requires the use of optical images of better quality. Considering 
the typical UV-optical colour magnitude relation of galaxies in 
the Virgo cluster (e.g. Boselli et al. 2014a), a cut at NUV - 
20 mag roughly corresponds to a cut in the SDSS g band of ~ 
19 mag for the blue star forming objects and g ~ 16 mag for 
the quiescent dwarf ellipticals. At this depth the SDSS is com¬ 
plete both in terms of total magnitude and surface brightness 
and can thus be used for this purpose. Indeed, out of the 33026 
GUViCS point-like sources with NUV < 20 mag, 32277 (98%) 
have an optical counterpart in the SDSS (the 749 objects without 
a SDSS counterpart are mainly ghosts of saturated stars, HII re¬ 
gions in resolved galaxies, or multiple systems unresolved in the 
UV images). About one third of these (7532) are photometrically 
identified as galaxies in the SDSS (SDSS_phot_ObjType=3). We 
have cross-matched these 7532 objects with the SDSS DR12 
spectroscopic catalogue, with NED, and with eBOSS through 
the DR12 SDSS finding chart tool. Spectra are available for 5620 
galaxies (75% of the sample). Their optical image on the SDSS 
DR 12 finding chart revealed, however, that a significant frac¬ 
tion of them are bright stars. Generally saturated, they are ex¬ 
tended on the SDSS image and thus missclassified as galaxies 
by the automatic SDSS pipeline. Among the bona fide galaxies 
of the UV-VPS catalogue of Voyer et al. (2014) with available 
redshift we identify 125 objects with vel < 3500 km s“^ Eive 

^ https://ned.ipac.caltech.edu/ 

^ https://www.sdss3.org/dr9/ 


of these, however, are left out of the sample because their NUV 
magnitude determined using our ad-hoc procedures for extended 
sources are below the limiting value of NUV - 20 mag. The re¬ 
sulting sample of Virgo cluster members determined using spec¬ 
troscopic data consists thus of 808 galaxies. 

3.2. Identification using optical NGVS scaling relations 

The identification of Virgo members among those objects with¬ 
out any spectroscopic data is critical in the determination of the 
luminosity function. The membership criteria adopted by the 
NGVS, fully described in Eerrarese et al. (2015), give the proba¬ 
bility that a galaxy belongs to Virgo according to its location in a 
multi-parameter space defined by a combination of galaxy struc¬ 
tural parameters, photometric redshifts (Raichoor et al. 2014), 
and an index measuring the strength of residual structures in im¬ 
ages created by subtracting from each galaxy a best fitting Sersic 
model. The exact combination of axes in this space was selected 
to allow for maximum separation of known (mostly spectroscop¬ 
ically confirmed) Virgo members and background sources, the 
latter identified in control fields located at three virial radii from 
the cluster centre. The NGVS covers 104 square degrees, and its 
entire footprint is mapped by GUViCS. We thus cross-matched 
the GUViCS UV-VPS catalogue with the NGVS catalogue of 
Virgo members identified using their photometric parameters 
and found 10 UV sources identified as Virgo members in the 
NGVS. We assume all these objects at a distance of 17 Mpc, 
the typical distance of Virgo and of its surrounding regions. The 
depth of the NGVS survey, which is 3 mag deeper than the SDSS 
in terms of surface brightness sensitivity, guarantees that all UV 
detected galaxies down to NUV < 20 mag have an optical coun¬ 
terpart. 

3.3. Identification using statistical estimates oustide the 
NGVS footprint 

The NGVS covers Virgo out to one virial radius, for a total of 
104 square degrees. To identify Virgo members in the GUViCS 
areas beyond the NGVS boundaries, and in the absence of spec¬ 
troscopic information, we must therefore resort to statistical ar¬ 
guments. In the GUViCS UV-VPS catalogue there are 1884 UV 
galaxies without any membership information. We drop those lo¬ 
cated within the NGVS footprint, reducing the sample of objects 
with no available membership information to 1155 sources. We 
then use three different statistical methods to estimate the num¬ 
ber of potential members among these galaxies. 

3.3.1. Statistical correction based on NGVS results 

If 10 galaxies with no redshift in the NGVS footprint (104 deg.^) 
are identified as Virgo members using criteria adopted by the 
NGVS, and the density of galaxies is conserved over the sur¬ 
veyed region, 20 extra objects are expected to be present in the 
remaining ~ 200 deg.^ covered by GUViCS. Since the NGVS 
observed the central region of the cluster, where the density is 
at its maximum, this number should be considered as an upper 
limit. 

3.3.2. Statistical correction based on redshift measurements 

Out of the 5620 galaxies in the VPS catalogue of Voyer et al. 
(2014) with spectra in the 300 deg.^ of the GUViCS survey only 
125 have vel < 3500 km s“*. If we assume that the same ratio of 
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members/background galaxies is conserved among those objects 
without any spectroscopic information (1884), we would thus 
expect 42 new members. This number drops to 30 if we limit this 
comparison in the region located outside the NGVS footprint 
corresponding to the cluster periphery (90/3452). 

3.3.3. Statistical correction based on GALEX-SDSS scaling 
relations 

An alternative method for identifying possible members has 
been first proposed by Boselli et al. (2011). For the same rea¬ 
sons given above, cluster galaxies populate different regions in 
the optical and UV surface brightness vs. observed magnitude 
relations and can thus be discriminated using their photometric 
parameters. Unfortunately the relations proposed by Boselli et 
al. (2011) are not sufficiently accurate to identify Virgo mem¬ 
bers among quiescent dwarf galaxies such as those located in 
the central 12 deg.^ analysed in that work. Here we propose an 
alternative technique based on three different scaling relations 
(S(r) vs. r mag; g - i vs. NUV - g, NUV - i vs. i mag, where 
S(r) is the mean surface brightness within the half-light radiufl0 
which we demonstrate to be fairly efficient at rejecting a signif¬ 
icant number of background objects also among star forming 
systems such as those detected in the outer regions of the cluster 
(see Fig.[T]i. Despite their relative low scatter, these relations are 
unfortunately not sufficiently accurate for an unambiguous iden¬ 
tification of cluster members and thus need to be combined with 
statistical tests. 

These scaling relations are determined using either the 
Extended Virgo Cluster Catalogue (EVCC, Kim et al. 2014) 
or the GUViCS catalogue of galaxies analysed in Boselli et al. 
(2014a). The EVCC catalogue has been chosen because it in¬ 
cludes galaxies at the same distance as Virgo (vel < 3000 km 
s '), it has both star forming and quiescent galaxies, and pro¬ 
vides structural parameters comparable to those listed in the 
GUViCS point-like catalogue of Voyer et al. (2014). The EVCC 
catalogue is used to determine the 2(r) vs. r apparent magni¬ 
tude relation. Unfortunately the EVCC catalogue does not in¬ 
clude UV magnitudes and thus cannot be used in the g - i vs. 
NUV - i colour-colour and the NUV - f vs. i relations. Eor these 
UV based relations we use as comparison sample the one com¬ 
posed of Virgo cluster galaxies analysed in Boselli et al. (2014a). 
Galaxies are separated according to their morphological type 
(early- from late-type) and atomic gas content. Eor consistency, 
the Virgo cluster sample is limited to NUV < 20 mag. 

Eigure[T] shows that; 

a) a number of sources identified by the standard SDSS pipeline 
or by our visual inspection of the SDSS finding charts as galax¬ 
ies are saturated stars. They have r-band surface brightnesses 
~ 5 mag brighter than those of typical galaxies of similar 
magnitude, but are rather comparable to those of stars. We 
identify as possible members in the 2(r) vs. r mag plot all 
galaxies satisfying the condition S(r) > 0.5xrH-12.4. We recall 
that this criterion is different and less strict than the one used in 
Boselli et al. (2014a) (the number of rejected objects is relatively 
small). This difference is due to the fact that in the inner 12 
deg.^ at faint magnitudes the population of Virgo members is 


The GUViCS catalogue lists the Petrosian half-light radius, while 
the EVCC the Kron half-light radius. We have checked that, for galaxies 
belonging to both samples, these radii are consistent. 

^ All values are corrected for Galactic extinction. These relations are 
also different from those adopted by Ferrarese et al. (2015) based only 
on optical data. 


strongly dominated by dwarf ellipticals with red colours. Thus 
the difference in S(r) with background star forming galaxies was 
more pronounced than in the present sample, which includes 
local star forming systems at the periphery of the cluster. 

b) all Virgo cluster members avoid the region above the g - i - 
0.6 X {NUV - g) - 0.2 and g-i^ 0.23 x {NUV - g) H- 0.45 
relations in the g - i vs. NUV - g colour-colour diagram. These 
regions are populated mainly by background sources. 

c) the sharp cut in the NUV - i vs. i colour-magnitude relation 
is due to the adopted cut in NUV. This relation, however, shows 
that all identified members have observed i magnitudes <18. 

These relations can be combined with a limit on the photo¬ 
metric redshift, accurately determined on these bright sources, to 
reject background galaxies. This is done by removing all objects 
with ZphotNN >0.1, where NN stands for the Neural Network 
(NN) photometric redshift definition (D’Abrusco et al. 2007). 
This cut in photometric redshift, which corresponds to the one 
adopted by Boselli et al. (2011) for the rejection of background 
galaxies in the inner 12 deg.^ of the cluster, has been chosen to 
maximise the number of background rejected objects and min¬ 
imise the exclusion of possible members. With this new selec¬ 
tion criteria, the sample of possible members outside the NGVS 
footprint drops to 264 objects. The number of spectroscopically 
confirmed Virgo galaxies satisfying the same photometric cri¬ 
teria is 75/90 (83.3%), while that of members over background 
objects 75/1636 (4.6%). We thus expect that 15 out of the 264 
galaxies satisfying the criteria are at the distance of Virgo. 

3.3.4. The applied statistical correction 

The three statistical methods described above give consistent re¬ 
sults. We thus include 15 extra objects to those already selected 
spectroscopically (see sect. 3.1) or using NGVS data (sect. 3.2) 
in the determination of the NUV luminosity function. The dis¬ 
tribution of the 264 galaxies selected according to Fig. [T]is used 
to make the appropriate statistical correction as a function of the 
NUV absolute magnitude in the determination of the NUV lumi¬ 
nosity function (see Fig. |2]l. These 15 galaxies are all assumed 
to be at a distance of 17 Mpc. Obviously the correction is more 
important at faint luminosities. Figure [1] also shows that these 
statistically added galaxies, all located in the outer regions of 
the cluster, have blue colours. 

The resulting sample is thus composed of 808 spectroscopi¬ 
cally confirmed sources, 10 objects identified at the distance of 
the Virgo cluster using NGVS based scaling relations, and 15 
extra objects added for statistical considerations, for a total of 
833 galaxies. This indicates that the determination of the NUV 
luminosity function is very robust against statistical uncertain¬ 
ties, making the Virgo cluster the ideal target for such studies in 
the nearby universe (e.g. Rines & Geller 2008). 

3.4. Possible biases 

The identification of galaxies through automatic pipelines is first 
done on surface brightness criteria. Low surface brightness, ex¬ 
tended systems are hardly detectable even if their integrated 
magnitude is above the limiting magnitude requested for the 
completeness of the survey. These systems might thus be poten¬ 
tially missed in the GUViCS survey. This effect, however, should 
be limited because 65% of the Virgo cluster region has been 
mapped at the depth of the GALEX MIS (~ 1500 sec), where 
the completeness limit of 22.7 mag (Morrissey et al. 2007) is 
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Fig. 1. The S(r) vs. r (upper row), g - i vs. NUV - g (middle row), and NUV - i vs. i (lower row) scaling relations used to 
identify Virgo cluster possible members. Left panels: scaling relations of spectroscopically confirmed Virgo members included in 
the extended source GUViCS catalogue. The black empty squares indicate the local EVCC galaxies. The red, cyan, and blue empty 
circles in the middle and lower panels indicate the early-type and late-type galaxies, these last selected according to their Hl-gas 
content (HI-def< 0.4, blue symbols; HI-def> 0.4, cyan symbols). The black solid lines indicates the limits in the different scaling 
relations used to identify possible Virgo cluster members. Central panels: relations drawn when all optically identified galaxies 
with NUV < 20 mag from the GUViCS VPS catalogue are included. Stars, plotted only in the central upper panel, are indicated 
with a brown open circle. VPS galaxies with vel < 3500 km s ' are indicated with magenta filled dots, spectroscopically confirmed 
background sources with filled yellow triangles, while potential members with photometric redshift <0.1 with green empty circles. 
Right panels: scaling relations of spectroscopically confirmed Virgo members (as in the left panels) with candidate members from 
the VPS catalogue selected using the criteria given in sect. 3.3.3 (green empty circles). All these relations are determined using UV 
data corrected for Galactic dust attenuation. 


significantly deeper than the one adopted in this work. At this 
depth low surface brightness features can be detected down to ~ 
28.5 mag arcsec^^ (e.g. Boselli et al. 2005b). Furthermore, Voyer 
et al. (2014) have systematically checked the UV images of all 
the optically selected VCC galaxies. The analysis presented in 
Boselli et al. (2011) revealed that most of the UV detected galax¬ 


ies used to construct the UV luminosity function down to NUV 
- 21 mag in the central 12 deg.^ of Virgo (1 mag below the UV 
limit adopted in this work) have a counterpart in the VCC cat¬ 
alogue. This central region, however, is dominated by quiescent 
galaxies with red colours, thus bright in the optical bands if de¬ 
tected in the NUV. 














6 


Boselli et al.: GUViCS: the UV luminosity function of the Virgo cluster 



Fig. 2. The NUV redshift completeness of the GUViCS sam¬ 
ple compared to that of Coma (red open squares; Cortese et al. 
2008), A1367 (green filled triangles; Cortese et al. 2005), and 
Shapley supercluster (orange filled squares; Haines et al. 2011). 
The vertical arrow shows the photometric completeness limit of 
the present survey (NUV - -11.45). 


None of the three extremely low surface brightness extended 
galaxies (2(g) ~ 28 mag arcsec^^) recently discovered by Mihos 
et al. (2015) in the core of the Virgo cluster has been detected by 
GALEX. Given their total magnitude (g > 17 mag), and the typ¬ 
ical colour of quiescent dwarf ellipticals (NUV-g > 3 mag), we 
would expect that their total NUV magnitude is > 20 mag, thus 
below the limiting magnitude adopted for the determination of 
the UV luminosity function. Less extreme objects such as those 
observed in Coma by van Dokkum et al. (2015) and Koda et 
al. (2015) are easily detected by the NGVS and are thus con¬ 
sidered in the NUV vs. optical cross-identification. The NGVS 
footprint largely overlaps with the VCC region, and has been ob¬ 
served in four photometric bands (ugiz) at a point-source depth 
of g ~ 25.9 mag and a surface brightness limit 2(g) ~ 29 mag 
arcsec^^ (Ferrarese et al. 2012). This limit in surface brightness 
is roughly 4 mag below the one reached by the VCC (Ferrarese et 
al. 2012). Among the UV detected galaxies with NUV < 20 mag, 
43 objects identified as Virgo members from spectroscopic mea¬ 
surements or using the NGVS scaling relations have an optical 
counterpart in the NGVS but not in the VCC. Outside the NGVS 
footprint optical data are available from the SDSS. Its sensitiv¬ 
ity in surface brightness is ~ 26 mag arcsec^^, thus ~ one mag 
better than the VCC. At the depth of the present survey (NUV 
< 20 mag) most of the UV detected sources have a counterpart 
in the SDSS (98%, see previous section), it is thus likely that 
the number of Virgo cluster galaxies without an optical counter¬ 
part missed in the present analysis is low. A rough estimate of 
this number can be derived thanks to the NGVS. If the density 
of galaxies is conserved all over the GALEX surveyed region, 
and if we assume that the sensitivity of the VCC is comparable 


to that of the SDSS, we expect to miss 86 galaxies (2 x 43). 
Once again this number should be considered as an upper limit 
because i) the density of galaxies significantly drops outside the 
central regions of the cluster, and ii) the SDSS is ~ 1 mag in 
surface brightness more sensitive than the VCC. Furthermore, 
this number is comparable to the number of galaxies statistically 
added in the previous section (15). We thus expect that the num¬ 
ber of faint galaxies missed for their low surface brightness in 
UV is small compared to the total number of galaxies. 

Another possible source of uncertainty is the presence of ex¬ 
tremely compact extragalactic sources erroneously identified as 
stars in the SDSS optical images whose typical resolution is of 
the order of ~ 1 arcsec. Typical examples are the ultra compact 
dwarf galaxies (UCDs), a population of extragalactic sources 
with characteristics in between those of compact elliptical galax¬ 
ies and globular clusters (Hilker et al. 1999; Drinkwater et al. 
2000; Zhang et al. 2015). These sources, however, have NUV 
apparent magnitudes below the NUV - 20 mag limit because of 
their relatively red colours (Boselli et al. 2011) and should not 
be counted in the determination of the NUV luminosity function 
of Virgo. Blue compact dwarf galaxies (BCD), characterised by 
an intense activity of star formation, have angular dimensions of 
>10 arcsec at the distance of Virgo and are thus fully resolved 
by the SDSS. They are thus included in the present analysis. 


4. The NUV luminosity function 

4.1. The NUV luminosity function of the whole cluster 

Figure |3] shows the NUV luminosity function of the whole 
cluster determined over the 300 deg.^ of the GUViCS survey 
The number of galaxies per bin of absolute NUV magnitude 
are given in Table [1] Given that the virial radius Ry of clus¬ 
ter A, the main substructure composing Virgo, is 5.38 degrees 
(McLaughlin 1999; Ferrarese et al. 2012, corresponding to 1.60 
Mpc at the distance of 17 Mpc), this luminosity function is rep¬ 
resentative of the cluster up to ~ 1.8 virial radii. For compari¬ 
son with other clusters we recall that the core radius of Virgo 
is of 0.45° (130 kpc; Boselli & Gavazzi 2006). We fit the data 
with a Schechter function with parameters given in Table |4] 
Following Ilbert et al. (2005) we use the MINUIT package of 
the CERN library (James & Roos 1995) to minimise the like¬ 
lihood (MIGRAD procedure), to obtain the non-parabolic error 
for each parameter (MINOS procedure), and the error contour a 
- M* (MNCONT procedure). The parameters of the fit are com¬ 
pared to those determined for the central 12 deg.^ of Virgo by 
Boselli et al. (2011), for other clusters and for the field in Fig.|4] 
The NUV luminosity function of Virgo is well fitted by a 
Schechter function with M* - -17.56 and a — -1.19. These val¬ 
ues are close to those determined in the central 12 deg.^ by 
Boselli et al. (2011), but are significantly different from those 
determined in Coma (Cortese et al. 2008, Hammer et al. 2012), 
A1367 (Cortese et al. 2005), in the Shapley supercluster (Haines 
et al. 2011), or in the field (Wyder et al. 2005). The slope of the 
faint end of the NUV luminosity function of Virgo (a = -L19)is 
significantly flatter than the one of Coma, A1367, or the Shapley 
supercluster (a -L5/-L6) as first noticed by Boselli et al. (2011). 
The slope is however comparable to the one determined in the 
field (a - -1.16; Wyder et al. 2005). M* is ~ 1 mag fainter than 
in Coma, in the Shapley supercluster, and in the field, and 2 mag 
fainter than in A1367. 


Corrected for Galactic attenuation. 
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Fig. 3. The NUV luminosity function of the Virgo cluster (filled 
dots). The best fit Schechter function (black solid line) is com¬ 
pared to the one determined for the central 12 deg.^. of the clus¬ 
ter by Boselli et al. (2011) at a limiting absolute magnitude of 
NUV < -10 (cyan dotted-dashed line) and NUV < -13 (cyan dot¬ 
ted line), and to that of the Coma cluster (red dotted line, from 
Cortese et al. 2008), A1367 (green long dashed line, Cortese 
et al. 2005), Shapley supercluster (orange dotted-dashed line; 
Haines et al. 2011) and the field (blue dashed line, Wyder et al. 
2005). The different luminosity functions have been normalised 
to ~ the same number of objects above the completeness limit of 
each cluster for a fair comparison. The vertical arrow shows the 
completeness limit of the present survey {NUV - -11.45). 


4.2. The luminosity function in the different cluster 
substructures 

The comparison of the properties of the luminosity function de¬ 
termined in different environments, from the cluster periphery to 
the densest regions in the core of the cluster, are useful to con¬ 
strain statistically the role played by the environment in shap¬ 
ing galaxy evolution. The determination of the statistical prop¬ 
erties of the different cluster substructures, accreted at different 
epochs, is important to understand the role of pre-processing in 
the transformation occurring in cluster galaxies. 

4.2.1. Cluster A vs. periphery 

Thanks to the large number of objects we can determine the 
NUV luminosity function separately for the core of the cluster 
and for its periphery. Consistently with Boselli et al. (2014a), we 
select galaxies within cluster A as those objects with vel < 3500 
km s"' and with an angular distance from M87 < 2.692°, corre¬ 
sponding to half its virial radius. The periphery is composed by 
all galaxies with vel < 3500 km s * and an angular distance from 
M87 > 6.1° (> 1.1 virial radii from the main body of the cluster). 
The NUV luminosity functions of cluster A and of the cluster 
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Fig. 4. Comparison of the parameters obtained after fitting a 
Schechter function to the NUV luminosity function of the Virgo 
cluster (black open circle), compared to those determined in 
Coma over 25 Mpc° ^ (red empty square, Cortese et al. 2008) 
or in its SW infalling region (red filled square. Hammer et al. 
2012), in A1367 (green filled triangle, Cortese et al. 2005), in the 
Shapley supercluster (orange filled square, Haines et al. 2011), 
and the field (blue asterisk, Wyder et al. 2005). Contours indicate 
1 <T probability distribution of the two correlated parameters. 


periphery, whose values are given in Table |2] are fitted with a 
Schechter function and compared in Fig. |5l while the fitted pa¬ 
rameters are given in Table |4] The fitted parameters for cluster A 
and for the periphery are compared to those determined for the 
whole cluster and for other clusters in Fig. |4] 

Figures |4] and |5] and Table |4] show that the NUV luminosity 
function, as in Coma (Cortese et al. 2008), does not significantly 
change from the core of cluster A to the cluster periphery outside 

> l-lRvir- 

4.2.2. The other cluster substructures 

The luminosity function of the main body of the cluster (cluster 
A) determined within half of the virial radius and that of 
the cluster periphery (> 1.1 Rvir) can also be compared to those 
determined for the other cluster substructures. This can be done 
only for cluster B and W cloud, which have a sufficient number 
of objects. To increase the statistics in the remaining substruc¬ 
tures, we combine their data according to the general properties 
of their members as determined in Boselli et al. (2014a). We 
thus combine cluster C and W’ cloud, which are dominated by 
quiescent objects, while the M cloud with the low velocity cloud 
(LVC) which, on the contrary, are mainly composed of star form¬ 
ing systems. The comparison of the different NUV luminosity 
functions is shown in the lower panel of Fig. |5] and in Fig. ??, 
while the fitted parameters are given in Table |4] Despite the poor 
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statistics, the comparison of the different luminosity functions 
suggests that these substructures have a flatter slope at the faint 
end than cluster A or the held. These trends, however, should be 
confirmed on stronger statistical basis. 

4.3. The luminosity function as a function of morphological 
type 

The sample can also be cut in different subclasses according to 
the morphological type of galaxies. This can be done either us¬ 
ing the accurate morphological classification given in the VCC 
(Binggeli et al. 1985) or in the EVCC (Kim et al. 2014), or 
using the colour of galaxies as done in Boselli et al. (2014a). 
Consistently with our previous work we select galaxies belong¬ 
ing to the red sequence, the green valley, and the blue cloud ac¬ 
cording to their location in the Mstar vs. NUV - i diagram us¬ 
ing the same selection criteria given in Boselli et al. (2014a). We 
prefer to select galaxies according to a quantitative colour-stellar 
mass relation rather than to a subjective morphology classifica¬ 
tion mainly because of the physical interpretation of this diagram 
in the framework of galaxy transformation in high-density envi¬ 
ronments (see Boselli et al. 2014a for details). 

Figure ?? shows the NUV luminosity function of the whole 
cluster separately for galaxies belonging to the red sequence, the 
green valley, and the blue cloud. The values are given in Table 
[1] while the parameters of the best fitted Schechter functions are 
given in Table |4] and compared in Fig. ??. Figures ?? and ??, 
and Tables [T] and 01 consistently indicate that at the bright end 
the NUV luminosity function of the Virgo cluster is dominated 
by late-type galaxies (the difference in M* between blue and red 
galaxies is of ~ 1 mag), while the increasing number of galaxies 
at the faint end is mainly due to the red quiescent galaxy pop¬ 
ulation, as already found in other clusters (Cortese et al. 2005, 
2008). Indeed the slope of the early-type red galaxies (a - -1.37) 
is significantly steeper than the one of the star forming systems 
within the blue cloud (a - -0.94). 

The comparison between the different morphological classes 
can also be done separately for galaxies located within cluster A 
and at the periphery of the cluster (Fig.|9j Table [TJ. The limited 
number of galaxies prevents an accurate determination of the 
parameters of the Schechter function for early-type galaxies in 
the periphery of the cluster. 

Despite a significant decrease of the total number of objects, 
the shape of the NUV luminosity function of red galaxies does 
not change significantly between cluster A and the periphery. 
A Kolmogorov-Smirnov test indicates that the probability that 
the cluster A and the cluster periphery early-type galaxy dis¬ 
tributions are driven by the same parent distribution is 11%. 
Important variations are instead observed in the NUV luminos¬ 
ity function of blue galaxies. Evident is the change in the slope 
at faint luminosities, from a = -1.10 in the outskirts of the clus¬ 
ter to a = -0.76 in the core of cluster A. In the cluster periphery 
the NUV luminosity function is dominated by blue galaxies at 
all luminosities (notice that the value of a = -1.10 determined 
for this category of objects is very close to the one derived in the 
held by Wyder et al. (2005)), while in the core of the cluster is 
dominated by the steep rise of red early-type dwarf galaxies be¬ 
low NUV > -14 mag (a - -1.25). This result does not necessary 
contradicts what found in the Coma cluster, i.e. that the shape 
of the NUV luminosity function of early- and late-type galax¬ 
ies does not change in the clustercentric distance range 0.5° < 
R < 2° (Cortese et al. 2008). Indeed, once considered in terms 
of virial radii, this range (0.3 < RjRvir < 1-2) does not match 


with the one considered in this analysis (cluster A; RjRvir < 0.5; 
cluster periphery;/?//?(/,> > 1.1). 

4.4. The star formation rate 

The UV emission of star forming systems is due to young and 
massive stars whose life on the main sequence is relatively short. 
If the star formation activity of galaxies is fairly constant over 
timescales as long as the typical age of the emitting stars (sta- 
tionarity conditions), UV luminosities corrected for dust attenu¬ 
ation can be converted into star formation rates (e.g. Kennicutt 
1998, Boselli et al. 2009; Boquien et al. 2014). The constant 
of transformation can be determined using different popula¬ 
tion synthesis models, and depends on the assumed IMF and 
metallicity (e.g. Boselli 2011). The NUV emission of late-type 
galaxies is principally due to stars with a typical age of ~ 150 
Myr (Boissier 2013). These timescales are comparable to the 
timescales for the observed variations in the star formation ac¬ 
tivity of cluster galaxies. Indeed, in these objects the interac¬ 
tion with the hostile cluster environment is able to remove, on 
relatively short timescales (~ 100-700 Myr, Boselli et al. 2006, 
2008a; Tonnesen & Bryan 2009), a large fraction of the atomic 
and molecular gas content (Fumagalli et al. 2009; Boselli et al. 
2014b) quenching the activity of star formation. This process is 
particularly efficient in dwarf systems, where the shallow poten¬ 
tial well of the galaxies cannot keep the gas reservoir anchored 
to the disc even in the central regions. The stationarity condi¬ 
tions required for transforming NUV luminosities into star for¬ 
mation rate are thus not necessarily satisfied. For comparison 
witfi otfier works and keeping in mind this caveat, however, we 
can try to convert the NUV luminosities corrected for dust atten¬ 
uation (see sect. 2) into SFR. Because of the reason mentioned 
above (quenching of the star formation activity of cluster galax¬ 
ies), using NUV data we expect to overestimate the present day 
star formation activity of cluster galaxies, in particular in low 
luminosity objects. 

Figure ?? shows the distribution of the star formation rate 
determined for late-type galaxies in the Virgo cluster region. 
This distribution has been determined for galaxies in the blue se¬ 
quence only (upper panel), or for all star forming galaxies (blue 
sequence and green valley; lower panel) separately for objects 
situated within cluster A and at the periphery of the cluster. We 
do not consider here red quiescent galaxies, where the UV emis¬ 
sion can be due to evolved stars (e.g. O’Connel 1999; Boselli 
et al. 2005b). The SFR luminosity distribution of both the blue 
sequence and of all the star forming galaxies at the periphery of 
the Virgo cluster is similar to the one determined for the field 
using SDSS and GAMA data (Gunawardfiana et al. 2013). It is, 
however, significantly steeper tfian tire one determined for galax¬ 
ies located witfiin cluster A, consistent witfi wfiat sfiown in Fig. 
|9] These distributions are fitted with a Schechter function whose 
parameters are given in Table 

5. Discussion 

The origins of the differences observed with other clusters and 
with the field (sect. 4.1) might be several, as already discussed 
in Boselli et al. (2011). They include: 

a) A different sampling in absolute magnitude for the different 
clusters and the field. The limiting magnitude in Coma, A1367, 
and Shapley is only NUV < -14.5 vs. ^ -11.5 in Virgo. If we 
limit the comparison with these clusters down to the limiting 
magnitude of NUV < -14.5 mag, however, we do not observe 
any significant change in the slope a of the Schechter function 
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Fig. 9. Upper panels: the NUV luminosity function of galaxies in the Virgo cluster periphery (left, filled dots) and in the core of 
cluster A (right, empty squares). Red symbols are for galaxies belonging to the red sequence, green symbols for galaxies in the 
green valley, and blue symbols for galaxies in the blue cloud. Lower panels: the NUV luminosity function of cluster A (empty 
squares) and cluster periphery (filled dots) galaxies in the blue (left) and red (right) sequences. The red and blue lines give the the 
best fit Schechter functions derived for the different samples. The different luminosity functions have been normalised to ~ the same 
number of objects. 


in Virgo. Given the degeneracy between the Schechter parame¬ 
ters a and M* it is difficult to establish whether the position of 
the knee of the Schechter function is significantly different in 
the various clusters and in the field or is just due to their steeper 
slope at the faint end. We recall, however, that clusters such as 
Coma and A1367 contain a number of gas-rich late-type sys¬ 
tems characterised by a strong activity of star formation proba¬ 
bly induced by their interaction with the surrounding intracluster 
medium (Iglesias-Paramo et al. 2004). These galaxies have been 
first discovered by Bothun & Dressier (1986) in the Coma clus¬ 
ter, and are typical of unrelaxed clusters such as A1367. Typical 


examples are the galaxy CGCG 97-73, 97-79, and 97-87 at the 
north west periphery of the cluster (Gavazzi et al. 2001). Similar 
objects are apparently lacking in Virgo as revealed by the com¬ 
plete Her narrow band imaging survey of the cluster (Koopmann 
et al. 2001; Boselli & Gavazzi 2002; Boselli et al. 2002; Gavazzi 
et al. 2006). 

b) Boselli et al. (2011) invoked a possible bias in the determi¬ 
nation of the UV luminosity function in the central 12 deg.^ of 
Virgo because of the lack of bright galaxies in their sample. The 
region analysed in this work is ~ 25 times larger and the number 
of objects has increased by a factor of ~ 8 with respect to that 
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work, making our determination of the NUV luminosity func¬ 
tion robust. 

c) Possible radial effects in the galaxy distribution. The Coma 
cluster NUV luminosity function of Cortese et al. (2008) has 
been determined over an area of 25 Mpc^ but avoiding the core 
of the cluster because of the presence of bright stars saturating 
the detector. On the contrary the NUV luminosity function of 
A1367 is limited to the central ~ 2 Mpc^, corresponding to ~ 
0.37 Ry. The properties of the fitted Schechter functions of these 
two clusters, however, are very similar, as depicted in Figs.[3]and 
m If limited to the infalling region in the south west of Coma, 
where UV data are available for galaxies down to NUV - -10.5, 
Hammer et al. (2012) determined a slope a - -1.39, slightly flat¬ 
ter than the one determined by Cortese et al. (2008) in the outer 
regions but still significantly steeper than the one determined in 
this work. Furthermore, we have shown that the shape of the fit¬ 
ted Schechter function does not change significantly from the 
core to the periphery of the cluster, corroborating the idea that 
radial effects are not important. 

d) The determination of the slope of the UV luminosity function 
at the faint end strongly depends on poorly constrained statisti¬ 
cal corrections in clusters other than Virgo (e.g. Rines & Geller 
2008; see Fig|2]l. At the limiting observed magnitude of 21 mag 
in the NUV used by Cortese et al. (2008) not all galaxies have 
an optical counterpart necessary for the star-galaxy identifica¬ 
tion on the SDSS (Voyer et al. 2014). In the Coma cluster and in 
the Shapley supercluster the completeness in redshift is present 
only down to NUV = -17, and drops to < 20% at the limiting 
absolute mag of NUV — -14 (Fig. |2]i. In Hammer et al. (2012) 
the UV selected galaxies (NUV < 24.5 mag) must have a SDSS 
optical counterpart with r < 21.2 mag. With this condition all 
star forming objects are automatically excluded by the sample. 
Although we cannot exclude real differences in the luminosity 
distribution of galaxies within these clusters, which are known to 
be characterised by different intrinsic properties (relaxation state 
and presence of substructures, spiral fraction, density of the intr¬ 
acluster medium and X-ray properties), it is conceivable that part 
of the observed differences results from the uncertain statistical 
corrections required in clusters more distant than Virgo. 

The NUV luminosity function determined in this work can 
also be compared to those obtained at other frequencies in the 
Virgo cluster. As summarised in Boselli & Gavazzi (2014), 
these are now available only in optical bands (Sandage et al. 
1985, photographic B-band; Trentham & Hodgkin 2002, B- 
band; Rines & Geller 2008, r-band; Lieder et al. (2012), V- and 
/-bands). We can compare our results to those obtained by Rines 
& Geller (2008) in the r-band within 1 Mpc from M87 since 
their sampled region fairly corresponds to the cluster A defini¬ 
tion adopted in this work (< 0.8 Mpc). Furthermore both works 
use membership criteria based on spectroscopic redshifts and are 
thus only marginally affected by statistical corrections. The fit 
with a Schechter function done by Rines & Geller (2008) gives 
= -21.32 and a - -1.28. Considering a typical NUV - r 
colour of ~ 5.5,4.5, and 3.5 mag for massive red sequence, green 
valley, and blue cloud galaxies as derived from the NUV - i 
colour-magnitude relation of Boselli et al. (2014a) and the typi¬ 
cal r - i colours of nearby galaxies (Fukujita et al. 1995), we can 
tentatively transform the NUV luminosity function of cluster A 
given in Table|4]in a r-band luminosity function. We thus expect 
that the derived r-band luminosity function is dominated by the 
early-type galaxy population with M* and a values close to those 
obtained by Rines & Geller (2008). We recall, however, that this 
is a very crude approximation since i) the colour of galaxies is 
luminosity-dependent and ii) the r- and NUV-bands are sensitive 


to very different stellar populations thus selection effects should 
be important. 

The variations of the NUV luminosity function with mor¬ 
phological type and cluster location can be analysed in the 
framework of galaxy evolution in high-density environments 
(see Boselli & Gavazzi 2006, 2014 for a review). The fact that 
the NUV luminosity function does not change significantly 
in terms of faint end slope going from the cluster core (< 0.5 
Ryir) to the outskirts (> 1.1 /?v,,>) and to the held, suggests that 
whatever perturbing mechanism is acting on low luminosity 
galaxies entering the cluster, this does not change their number 
per unit NUV luminosity. The interpretation of this observa¬ 
tional result, however, is not straightforward since the different 
environmental processes can modify in a complex way both the 
number of objects and the luminosity of the perturbed galaxies. 
As a first indication we would be tempted to exclude for Virgo 
the combined effect of dwarf tidal disruption in the core of the 
cluster and galaxy harassment invoked by Popesso et al. (2006), 
Barkhouse et al. (2007), and de Filippis et al. (2011) to explain 
the observed flattening in the core and the steepening in the 
periphery in the optical luminosity function of nearby clusters. 
Indeed, we do not observe in Virgo any radial variations of the 
luminosity function both in the NUV and in the r-band (Rines & 
Geller 2008). On the contrary, our results seem consistent with 
a mild transformation of star forming systems into quiescent 
dwarf elliptical galaxies after a rapid ram pressure stripping 
event as first proposed by Boselli et al. (2008a,b). Because of 
their shallow potential well low luminosity star forming galax¬ 
ies easily lose their gas because of ram pressure stripping and 
quench their activity of star formation on timescales of < 1 Gyr. 
This happens at the periphery of the cluster once they encounter 
the intergalactic medium trapped within the potential of the 
cluster. Because this timescale is relatively short compared 
to the crossing time of the cluster (~ 1.7 Gyr), these newly 
accreted galaxies are already red when they reach the cluster 
core (Boselli et al. 2014a). This effect is consistent with the 
observed reversed shape of the NUV luminosity function of red 
and blue galaxies in the cluster core and at the periphery (Fig|9]l. 
It is also consistent with the spectrophotometric, structural 
(Boselli et al. 2008a,b), and kinematic (Toloba et al. 2009, 2011, 
2012, 2015) properties of Virgo cluster dwarf elliptical galaxies, 
as extensively reviewed in Boselli & Gavazzi (2014). 

The perturbations induced by the cluster environment on the 
star formation activity of galaxies are obviously more efficient 
in low-mass systems than in massive spirals. In massive objects 
the steep potential well can retain some gas in the inner regions 
producing radially truncated discs as in the case of NGC 4569 
(Boselli et al. 2006). The timescale necessary for the total abla¬ 
tion of the atomic and molecular gas here is longer than in dwarfs 
as indicated both by observations and hydrodynamic simulations 
(Roediger & Bruggen 2007, Tonnesen & Bryan 2009, Hughes & 
Cortese 2009, Cortese & Hughes 2009, Fumagalli et al. 2009, 
Gavazzi et al. 2013a,b, Boselli et al. 2014a,b). The observed 
increasing difference with decreasing star formation rate in the 
S FR luminosity function of the Virgo cluster shown in Fig. ?? 
is a further statistical evidence of this phenomenon. 

6. Conclusion 

We have determined the properties of the NUV luminosity func¬ 
tion of the Virgo cluster over ~ 300 deg.^ (~ up to 1.8 virial radii) 
thanks to the GUViCS survey using a complete sample of ~ 833 
galaxies, out of which 808 identified as Virgo members using 
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spectroscopic data and 10 with optical scaling relations, down to 
the 100% completeness limit of NUV - -11.5. The results of 
this analysis can be summarised as follows: 

1) The NUV luminosity function can be represented by a 
Schechter function with M* =-17.56 and a =-1.19. M* is ~ 1-2 
mag fainter and the slope signihcantly flatter than the one found 
in other clusters. M* is only ~ 0.8 mag fainter and a compara¬ 
ble to those determined in the field. The difference with other 
clusters is probably due to the quite uncertain statistical correc¬ 
tion necessary to identify members in clusters at distances higher 
than Virgo. 

2) The NUV luminosity function of galaxies of different mor¬ 
phological type as determined according to their position within 
the NUV—i vs. Mstar relation is signihcantly different. Late-type 
galaxies belonging to the blue sequence dominate the bright end 
of the luminosity function, while their faint end is very Hat a - 
-0.94). On the contrary, early-type galaxies belonging to the red 
sequence dominate the faint end {a - -1.37) of the total NUV 
luminosity function below NUV ^ -13 mag. 

3) The overall properties of the NUV luminosity function de¬ 
termined within 0.5 virial radii of cluster A, the main structure 
of Virgo, and in the cluster periphery, outside ~ 1.1 virial radii, 
are comparable. In the outskirts of the cluster the NUV luminos¬ 
ity function is comparable to the one determined in the Held by 
Wyder et al. (2005). Here it is dominated down to NUV ^ -11.5 
by blue galaxies, with a mild rise at faint luminosities (a=-1.10). 
In the cluster core red quiescent galaxies dominate for NUV > 
-15 mag (a - -1.25), while the number of star forming systems 
signihcantly drops (a - -0.76). These observed properties of the 
NUV luminosity function are consistent with a rapid transfor¬ 
mation of star forming galaxies recently entered in the cluster 
environment after a ram-pressure stripping event. 

Acknowledgements. We thank the anonymous referee for constructive com¬ 
ments that helped improving the quality of the manuscript. This reseai'ch has 
been financed by the French ANR grant VIRAGE and the French national pro¬ 
gram PNCG. We wish to thank the GALEX Time Allocation Committee for the 
generous allocation of time devoted to this project. LC acknowledges financial 
support from the Australian Research Council (DP130100664). This research 
has made use of the NASA/IPAC Extragalactic Database (NED) which is oper¬ 
ated by the Jet Propulsion Laboratory, California Institute of Technology, un¬ 
der contract with the National Aeronautics and Space Administration and of 
the GOLDMine database I http://goldmine.mib.infn.it/ 1 (Gavazzi et al. 2003). 
GALEX (Galaxy Evolution Explorer) is a NASA Small Explorer, launched in 
2003 April. We gratefully acknowledge NASA support for construction, oper¬ 
ation, and science analysis for the GALEX mission, developed in cooperation 
with the Centre National d’Etudes Spatiales of France and the Korean Ministry of 
Science and Technology. Funding for the SDSS and SDSS-II has been provided 
by the Alfred P. Sloan Foundation, the Participating Institutions, the National 
Science Foundation, the U.S. Department of Energy, the National Aeronautics 
and Space Administration, the Japanese Monbukagakusho, the Max Planck 
Society, and the Higher Education Funding Council for England. The SDSS 
Web Site is http://www.sdss.org/ The SDSS is managed by the Astrophysical 
Research Consortium for the Participating Institutions. The Participating 
Institutions are the American Museum of Natural History, Astrophysical Institute 
Potsdam, University of Basel, University of Cambridge, Case Western Reserve 
University, University of Chicago, Drexel University, Fermilab, the Institute 
for Advanced Study, the Japan Participation Group, Johns Hopkins University, 
the Joint Institute for Nuclear Astrophysics, the Kavli Institute for Particle 
Astrophysics and Cosmology, the Korean Scientist Group, the Chinese Academy 
of Sciences (LAMOST), Los Alamos National Laboratory, the Max-Planck- 
Institute for Astronomy (MPIA), the Max-Planck-Institute for Astrophysics 
(MPA), New Mexico State University, Ohio State University, University of 
Pittsburgh, University of Portsmouth, Princeton University, the United States 
Naval Obseiwatory, and the University of Washington. 

References 

Abazajian, K. N., Adelman-McCarthy, J. K., Agiieros, M. A., et al. 2009, ApJS, 
182, 543 


Barkhouse, W. A., Yee, H. K. C., & Lopez-Cruz, O. 2007, ApJ, 671, 1471 
Binggeli, B., Sandage, A., & Tammann, G. A. 1985, AJ, 90, 1681 
Binggeli, B., Sandage, A., & Tammann, G. A. 1988, ARA&A, 26, 509 
Boissier, S. 2013, Planets, Stars and Stellar Systems. Volume 6: Extragalactic 
Astronomy and Cosmology, 141 

Boissier, S., Boselli, A., Due, P.-A., et al. 2012, A&A, 545, A142 
Boissier, S., Boselli, A., Voyer, E., et al., 2015, arXiv: 1504.06111 
Boquien, M., Buat, V., & PeiTet, V. 2014, A&A, 571, AA72 
Boselli, A. 2011, A Panchromatic View of Galaxies, by Alessandro Boselli. - 
Practical Approach Book - ISBN-10: 3-527-40991-2. ISBN-13: 978-3-527- 
40991-4 - Wiley-VCH, Berlin 2011. XVI, 324pp, Hardcover, 

Boselli, A., & Gavazzi, G. 2002, A&A, 386, 124 
Boselli, A., & Gavazzi, G. 2006, PASP, 118, 517 
Boselli, A., & Gavazzi, G. 2014, A&A Rev., 22, 74 

Boselli, A., Iglesias-Paramo, J., Vilchez, J. M., & Gavazzi, G. 2002, A&A, 386, 
134 

Boselli, A., Cortese, L., Deharveng, J. M., et al. 2005a, ApJ, 629, L29 

Boselli, A., Boissier, S., Cortese, L., et al. 2005b, ApJ, 623, L13 

Boselli, A., Boissier, S., Cortese, L., et al. 2006, ApJ, 651, 811 

Boselli, A., Boissier, S., Cortese, L., & Gavazzi, G. 2008a, ApJ, 674, 742 

Boselli, A., Boissier, S., Cortese, L., & Gavazzi, G. 2008b, A&A, 489, 1015 

Boselli, A., Boissier, S., Cortese, L., et al. 2009, ApJ, 706, 1527 

Boselli, A., Boissier, S., Heinis, S., et al. 2011, A&A, 528, A107 

Boselli, A., Voyer, E., Boissier, S., et al. 2014a, A&A, 570, AA69 

Boselli, A., Cortese, L., Boquien, M., et al. 2014b, A&A, 564, A67 

Bothun, G. D., & Dressier, A. 1986, ApJ, 301, 57 

Byrd, G., & Valtonen, M. 1990, ApJ, 350, 89 

Cortese, L., & Hughes, T. M. 2009, MNRAS, 400, 1225 

Cortese, L., Boselli, A., Gavazzi, G., et al. 2005, ApJ, 623, L17 

Cortese, L., Gavazzi, G., & Boselli, A. 2008, MNRAS, 390, 1282 

Cowie, L. L., & Songaila, A. 1977, Nature, 266, 501 

D’Abrusco, R., Staiano, A., Longo, G., et al. 2007, ApJ, 663, 752 

Davies, J. I., Baes, M., Bendo, G. J., et al. 2010, A&A, 518, L48 

de Filippis, E., Paolillo, M., Longo, G., et al. 2011, MNRAS, 414, 2771 

Dressier, A. 1980, ApJ, 236, 351 

Dressier, A. 2004, Clusters of Galaxies: Probes of Cosmological Structure and 
Galaxy Evolution, 206 

Dressier, A., Oemler, A., Jr., Couch, W. J., et al. 1997, ApJ, 490, 577 
Drinkwater, M. J., Jones, J. B., Gregg, M. D., & Phillipps, S. 2000, PASA, 17, 
227 

Ferrarese, L., Cote, R, Cuillandre, J.-C., et al. 2012, ApJS, 200, 4 
Ferrarese, L., Cote, R, MacArthur, L., et al., 2015, submitted to ApJS 
Fitzpatrick, E. L., & Massa, D. 2007, ApJ, 663, 320 
Fukugita, M., Shimasaku, K., & Ichikawa, T. 1995, PASP, 107, 945 
Fumagalli, M., Krumholz, M. R., Prochaska, J. X., Gavazzi, G., & Boselli, A. 
2009, ApJ, 697, 1811 

Gavazzi, G., Boselli, A., Scodeggio, M., Pierini, D., & Belsole, E. 1999, 
MNRAS, 304, 595 

Gavazzi, G., Boselli, A., Mayer, L., et al. 2001, ApJ, 563, L23 
Gavazzi, G., Boselli, A., Donati, A., Franzetti, R, & Scodeggio, M. 2003, A&A, 
400,451 

Gavazzi, G., Boselli, A., Cortese, L., et al. 2006, A&A, 446, 839 
Gavazzi, G., Fumagalli, M., Fossati, M., et al. 2013a, A&A, 553, A89 
Gavazzi, G., Savorgnan, G., Fossati, M., et al. 2013b, A&A, 553, A90 
Giovanelli, R., Haynes, M. R, Kent, B. R., et al. 2005, AJ, 130, 2598 
Gunawardhana, M. L. R, Hopkins, A. M., Bland-Hawthorn, J., et al. 2013, 
MNRAS, 433, 2764 

Gunn, J. E., & Gott, J. R., Ill 1972, ApJ, 176, 1 

Haines, C. R, Busarello, G., Merluzzi, R, et al. 2011, MNRAS, 412, 127 
Hammer, D. M., Homschemeier, A. E., Salim, S., et al. 2012, ApJ, 745, 177 
Hao, C.-N., Kennicutt, R. C., Johnson, B. D., et al. 2011, ApJ, 741, 124 
Hilker, M., Infante, L., Vieira, G., Kissler-Patig, M., & Richtler, T. 1999, A&AS, 
134,75 

Hughes, T. M., & Cortese, L. 2009, MNRAS, 396, L41 
Kennicutt, R. C., Jr. 1998, ARA&A, 36, 189 
Kim, S., Rey, S.-C., Jerjen, H., et al. 2014, ApJS, 215, 22 
Koda, J., Yagi, M., Yamanoi, H., & Komiyama, Y. 2015, ApJ, 807, L2 
Koopmann, R. A., Kenney, J. D. R, & Young, J. 2001, ApJS, 135, 125 
Iglesias-Paramo, J., Boselli, A., Gavazzi, G., & Zaccardo, A. 2004, A&A, 421, 
887 

Ilbert, O., Tresse, L., Zucca, E., et al. 2005, A&A, 439, 863 
James, R, Roos, M., 1995, MINUIT Function Minimization and Error Analysis, 
Version 95.03, CERN Program Library D506 
Larson, R. B., Tinsley, B. M., & Caldwell, C. N. 1980, ApJ, 237, 692 
Lieder, S., Lisker, T., Hilker, M., Misgeld, I., & Durrell, P. 2012, A&A, 538, A69 
McLaughlin, D. E. 1999, ApJ, 512, L9 
Mei, S., Blakeslee, J. R, Cote, R, et al. 2007, ApJ, 655, 144 
Merritt, D. 1983, ApJ, 264, 24 



12 


Boselli et al.: GUViCS: the UV luminosity function of the Virgo cluster 


Mihos, J. C., Durrell, P. R., Ferrarese, L., et al. 2015, ApJ, 809, L21 
Moore, B., Lake, G., & Katz, N. 1998, ApJ, 495, 139 
Monissey, P., Conrow, T., Barlow, T. A., et al. 2007, ApJS, 173, 682 
Nulsen, P. E. J. 1982, MNRAS, 198, 1007 
O’Connell, R. W. 1999, ARA&A, 37, 603 

Popesso, R, Biviano, A., Bohiinger, H., & Romaniello, M. 2006, A&A, 445, 29 

Raichoor, A., Mei, S., Erben, T., et al. 2014, ApJ, 797, 102 

Rines, K., & Geller, M. J. 2008, AJ, 135, 1837 

Roediger, E., & Briiggen, M. 2007, MNRAS, 380, 1399 

Sandage, A., Binggeli, B., & Tammann, G. A. 1985, AJ, 90, 1759 

Schlegel, D. J., Finkbeiner, D. R, & Davis, M. 1998, ApJ, 500, 525 

Toloba, E., Boselli, A., Gorgas, J., et al. 2009, ApJ, 707, L17 

Toloba, E., Boselli, A., Cenarro, A. J., et al. 2011, A&A, 526, A114 

Toloba, E., Boselli, A., Peletier, R. R, et al. 2012, A&A, 548, A78 

Toloba, E., Guhathakurta, R, Boselli, A., et al. 2015, ApJ, 799, 172 

Tonnesen, S., & Bryan, G. L. 2009, ApJ, 694, 789 

Trentham, N., & Hodgkin, S. 2002, MNRAS, 333, 423 

van Dokkum, P. G., Romanowsky, A. J., Abraham, R., et al. 2015, ApJ, 804, L26 
Voyer, E. N., Boselli, A., Boissier, S., et al. 2014, A&A, 569, AA124 
Whitmore, B. C., Gilmore, D. M., & Jones, C. 1993, ApJ, 407, 489 
Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868 
Wyder, T. K., Treyer, M. A., Milliard, B., et al. 2005, ApJ, 619, L15 
Zhang, H.-X., Peng, E. W., Cote, R, et al. 2015, ApJ, 802, 30 


1000 


fllOO 


10 - 


1 - 


>©' 


^100 


10 - 


1 - 


e- 


0.1 


T—I—I—I—I—I—I—I—I—I—I—I—I—I—I—I—I—I—I—I—rg 


r 


□ Cluster A 
A VCC/NGVS 


I.- ■ i' 








.1 yf 


y r 


!■ -I I I I I I I I I I I M l I I . 


H Cluster B 
^ W cloud 
☆ M + LVC clouds 
O C + W clouds 


32 23 17Mpc 


•r-f 




-L.— V- j: - 


AL 


/ 

i f 

I :• I 


'f/ 


n .;/■ 


/■' 

//■ 

// 


M <) 


T 


J_I_I _I_ I_I_L 




-20 -18 -16 -14 -12 

NUV (AB mag) 


Fig. 5. Upper panel; the NUV luminosity function of Virgo clus¬ 
ter A (magenta empty squares) and of the cluster periphery (cyan 
filled dots) are compared to the one determined over the 104° ^ 
sampled by the VCC and the NGVS survey (black filled trian¬ 
gles). The best fit Schechter functions are shown by the magenta 
dotted, cyan long dashed, and black dotted lines, respectively. 
The vertical arrow shows the limit of completeness of the lumi¬ 
nosity function. Lower panel; the same cluster A and periphery 
luminosity functions are compared to those determined for clus¬ 
ter B (grey empty stars), W cloud (brown filled triangles), M 
plus Low Velocity Clouds (dark green empty stars), and cluster 
C plus W’ cloud (orange empty hexagons). The vertical arrows 
show the different limits of completeness for the different clus¬ 
ter substructures; 17 Mpc for cluster A, cluster C, LVC, and the 
cluster periphery, 23 Mpc for cluster B and W’ cloud, 32 Mpc 
for W and M clouds. 
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Table 1. The NUV luminosity function. 


NUV 

All 

Blue 

Green 

Red 

-20 

- 

- 

- 

- 

-19 

5 

5 

- 

- 

-18 

14 

14 

- 

- 

-17 

59 

48 

2 

9 

-16 

68 

45 

13 

10 

-15 

127 

69 

18 

40 

-14 

156 

72 

36 

48 

-13 

196 

70 

47 

79 

-12 

204 

49 

53 

102 


Note: number of objects 

Table 2. The NUV luminosity function in the different cluster substructures. 


NUV 

Cluster A Cluster B 

W cloud 

C-rW’ 

M + LVC 

Periphery 

-20 

- 

- 

- 

- 

- 

- 

-19 

- 

1 

- 

- 

- 

1 

-18 

2 

- 

- 

- 

- 

4 

-17 

13 

5 

1 

3 

2 

12 

-16 

9 

10 

- 

4 

4 

21 

-15 

25 

6 

7 

7 

6 

36 

-14 

26 

10 

12 

4 

9 

46 

-13 

34 

13 

19 

8 

14 

62 

-12 

38 

11 

14 

6 

6 

63 

Note: number of objects 






Table 3. 

The NUV luminosity function; 

cluster A vs. periphery 




NUV Cluster A red 

Cluster A green 

Cluster A blue 

Periphery red 

Periphery green 

Periphery blue 

-20 

- 

- 

- 

- 

- 

- 

-19 

- 

- 

- 

- 

- 

1 

-18 

- 

- 

2 

- 

- 

4 

-17 

3 

2 

8 

2 

- 

10 

-16 

2 

1 

6 

- 

5 

16 

-15 

12 

6 

7 

6 

1 

29 

-14 

12 

6 

8 

6 

8 

32 

-13 

22 

7 

5 

9 

12 

41 

-12 

26 

8 

4 

17 

14 

32 


Note: number of objects 
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Table 4. The best fitting parameters for the NUV luminosity function. 


Sample 

a 

M* 

0. 

Full sample 

_1 1Q+U.UJ 

^•^^-0.03 

-17 56‘^o-23 

i/.JU-0 24 

88.40!“;” 

Red sequence 

1 '17+0-12 
-0,10 

-16 57‘^o.70 
-0.98 

94 93+20.59 
^^•^'^-13.35 

Green valley 

171 +0.10 
^•^^-0,10 

-16 11+0-40 

27.74!‘'4f 

Blue cloud 

-0,05 

-17 51‘^0-25 
^ ^•‘^^-0.25 

85.24![||“ 

Late-type 

^•^'-0,04 

-17 54+0.24 
i / .JH--0 25 

85.58!||” 

Cluster A 

1 77+0.08 
^•^^-0,08 

-17 49+0-50 
^ ^-^^-0.58 

14 14+6.37 
-5.19 

Cluster periphery 

1 77+0.07 
^•^^-0.07 


25 51+10-45 
z.j.j1-8_03 

Cluster A red sequence 

1 7C+O.24 
^•^^-0.28 

-15.70!“;” 

12.92!”;” 

Cluster A green valley 


-15.23!“;” 

10.8l!“;g 

Cluster A blue could 

-0 76+^-^^ 

-17.29!“;5“ 

14.62!““ 

Cluster A late-type 

-0 99+^1‘ 
^•^^-0,11 

-17 41+0-50 
^ -0.62 


Cluster periphery red sequence 

^•^^-0,16 

-23.97!],”^ 

0 12+4-46 
u.iZ,-oi2 

Cluster periphery green valley 


-13.64!“;58 

70 07+15.14 
Zo.oZ_j.^ 29 

Cluster periphery blue cloud 

-1 10+0-09 
^•^^-0,08 

-17.32!“;“ 

24.99!‘‘2“ 

Cluster periphery late-type 

1 1 c+0.08 
^•^^-0.07 

1 7 77+0.50 
^ ^-^'-0.53 

27.13!>”i 

Cluster B 

-1 i2‘^o.ii 

-18.41!“:“^ 

6.54”“ 

W cloud 

^•^^-0,17 

-15.94”;” 

1 9 77+36.81 
' -6.37 

Cluster C + W’ cloud 

i.oj-o 19 


C 70+6.10 
^•^°-4.91 

M cloud + LVC 

-0 93+0-24 

-16.22!“;“^ 

11.61”“ 

NGVS footprint'' 

1 1 q+0.04 
^•^°-0,04 

-17 

i / .jU-0 27 

58-78!”:°^ 

NGVS footprint" red sequence 

-1.36t|^ 

-16 31‘^^-”^'^ 

21 57‘^^o.9i 
-12.58 

NGVS footprint" late-type 

-0.98!“;“ 

-17 36‘^^-^^ 
i/.j»U-0 28 

64 92+15-60 

-12.50 


Note: <3=1 virial radius from M87 and M49. 


Table 5. The SFR luminosity function. 


log SFR 

Cluster A blue-l-green 

Periphery blue-)-green 

Cluster A blue 

Periphery blue 

-2.5 

12 

68 

4 

51 

-2.0 

17 

52 

10 

41 

-1.5 

15 

44 

7 

39 

-1.0 

6 

24 

1 

21 

-0.5 

10 

17 

9 

13 

0 

7 

11 

6 

9 

0.5 

4 

4 

4 

4 

1.0 

0 

1 

0 

1 


Note: number of objects 

Table 6. The best fitting parameters for the SFR luminosity function. 


Sample 

a 

logSFR* 

tp. 

Cluster periphery, blue cloud 

-1 26 ^ 0-05 

^•^^- 0.05 

0 54 + 0.2 / 

^•'^^- 0.31 

Q qt + 4.66 
®-®^- 3.09 

Cluster periphery, all star forming 

1 7 Q + 0.05 
^•^°- 0.04 

- 0.27 

10.26!^“ 

Cluster A, all star forming 

111 + 0.08 
- 0.08 

0.83!“:” 

6.58”““ 
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